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ABSTRACT 

In  this  paper,  we  investigate  the  kinematics  and  dynamics  of  floating , 
planar  four-bar  linkages.  The  geometry  of  configuration  space  is  analyzed 
through  the  classical  theory  of  mechanisms  due  to  Grashof.  The  techniques 
of  symplectic  and  Poisson  reduction  are  used  to  understand  the  dynamics 
of  the  system.  Bifurcations  of  relative  equilibria  for  linkages  admitting 
symmetric  shapes  are  studied  using  the  techniques  of  singularity  theory.  The 
problem  of  reconstruction  of  the  full  dynamics  and  its  relation  to  geometric 
phases  is  discussed  through  some  examples.  This  research  reveals  that  a 
coupled  mechanical  system  with  kinematic  loops  possesses  richer  and  more 
complicated  dynamical  aspects  in  comparison  with  systems  which  have  the 
same  number  of  degrees  of  freedom,  but  no  kinematic  loops. 


*  This  work  was  supported  in  part  by  the  AFOSR  University  Research  Initiative  Pro¬ 
gram  under  grant  AFOSR-90-0105  and  by  the  National  Science  Foundation’s  Engineering 
Research  Centers  Program:  NSFD  CDR  SS03012. 
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1.  INTRODUCTION 

There  has  been  significant  progress  in  our  understanding  of  the  hamiltonian  structure 
of  serial-link  (or  open  chain)  multibody  systems  (Baillieul  [1988];  Baillieul  Sz  Levi  [1987]; 
Grossman,  Krishnaprasad  &  Marsden  [1988];  Krishnaprasad  Sz  Marsden  [1987];  Marsden, 
Krishnaprasad  Sz  Simo  [1988];  Oh  et  al.  [1989];  Posbergh  [1988];  Posbergh,  Krishnaprasad 
8z  Marsden  [1987];  Simo,  Marsden  Sz  Krishnaprasad  [1988];  Simo,  Posbergh  Sz  Marsden 
[1989];  Sreenath  [1987];  Sreenath  et  al.  [1988];  Wang  Sz  Krishnaprasad  [1989]).  The 
application  of  geometric  methods,  symmetry  principles  and  reduction  has  led  to  deeper 
knowledge  of  the  dynamics  of  model  problems.  This  insight  has  also  been  helpful  in 
developing  appropriate  control-theoretic  tools  (Bloch,  Reyhanoglu  Sz  McClamroch  [1991]; 
Krishnaprasad  Sz  Yang  [1991];  Sreenath  [1990]).  The  primary  source  of  motivation  for 
such  problems  has  been  in  aerospace  engineering  where  imaginative  designs  of  multibody 
spacecraft  have  been  proposed  and  on  occasion  realized  (Bianchi  Sz  Schielen  [1985];  CIME 
[1972];  Wittenburg  [1977]). 

On  the  other  hand,  problems  in  aerospace  engineering  also'  suggest  that  multibody 
systems  with  kinematic  loops  are  of  practical  importance.  An  important  example  may 
be  the  parallel  linkage  based  robot  manipulators  which  have  been  contemplated  for  space 
applications.  However,  as  we  shall  see,  the  presence  of  loop  closure  constraints  implies 
that  the  knowledge  of  hamiltonian  structure  and  phase  portraits  for  open  chain  multibody 
systems  cannot  be  applied  directly  to  systems  with  closed  loops.  For  instance,  for 
multibody  systems  with  closed  loops,  the  topological  description  of  configuration  space 
strongly  depends  on  the  choice  of  kinematic  parameters.  This  is  not  true  for  systems  with 
open  chains. 

In  this  paper,  we  study  the  kinematics  and  dynamics  of  the  simplest  mechanical 
system  with  a  kinematic  loop  —  the  planar,  floating  four-bar  linkage.  By  floating  we 
mean  that  no  link  is  fixed  in  space.  It  is  simplest  in  that  it  has  the  fewest  degrees  of 
freedom  (after  symmetry  reduction)  among  all  kinds  of  mechanisms  with  closed  loops. 
In  the  theory  of  machines,  the  four-bar  linkage  (with  one  link,  the  frame,  fixed)  plays 
an  extremely  important  role  (Hunt  [197S]).  It  may  be  directly  used  in  many  mechanical 
systems.  More  complex  mechanisms  are  often  synthesized  using  one  four-bar  linkage  to 
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drive  another.  It  is  of  further  interest  that  a  four-bar  linkage  can  generate  the  wide  variety 
of  motions  which  are  represented  by  coupler  curves  (Hrones  &  Nelson  [1951]).  For  floating 
four-bar  linkages,  one  expects  similar  properties  to  hold.  As  we  shall  show  in  this  paper, 
floating  four-bar  linkages  do  possess  very  interesting  features  from  both  kinematic  and 
dynamic  points  of  view. 

The  outline  of  this  paper  is  as  follows.  After  stating  the  basic  notations  for  this  paper 
in  section  2,  we  give  a  geometric  description  of  the  configuration  space  in  section  3.  In 
section  4,  we  give  an  explicit  expression  for  the  kinetic  energy  which  is  also  the  lagrangian 
since  in  this  paper  we  assume  that  no  potential  energy  is  involved.  From  this  section 
we  see  that  explicit  computation  of  constrained  dynamics  is  difficult.  Yet  in  the  present 
setting,  using  geometric  techniques,  one  can  infer  qualitative  properties  without  recourse 
to  explicit  analytic  representation  of  the  constrained  dynamics.  This  we  do  in  section  5, 
wherein  we  explore  symmetry  properties,  hamiltonian  structure  and  reduction  of  four  bar 
linkage  dynamics.  In  section  6  we  use  a  theorem  of  Smale  to  compute  relative  equilibria 
for  the  dynamics  of  a  four-bar  linkage.  Further,  using  techniques  from  singularity  theory, 
we  study  the  local  bifurcations  of  relative  equilibria  for  linkages  which  admit  symmetric 
shapes.  In  section  7,  the  phase  reconstruction  problem  for  a  four-bar  linkage  is  investigated. 
It  is  interesting  to  note  that  it  is  possible  to  generate  either  positive  or  negative  phase  shift 
by  the  motion  in  shape  space.  This  phenomenon  does  not  arise  in  a  planar  two-body 
system  which  also  has  two  degrees  of  freedom. 

Preliminary  versions  of  this  paper  were  presented  at  the  IEEE  Conference  on  Decision 
and  Control  (Yang  &  Krishnaprasad  [1989];  Yang  &  Krishnaprasad  [1990]). 

We  are  pleased  to  acknowledge  useful  discussions  on  the  topics  of  this  paper  with 
Jerrold  Marsden  and  Tim  Healey. 
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2.  NOTATIONS  AND  GEOMETRIC  CONSTRAINTS 

The  structure  of  a  closed  floating  four-bar  linkage  is  represented  in  Figure  2.1.  What 
we  mean  by  bar  in  this  paper  is  a  planar  rigid  body,  on  which  the  center  of  mass  and  pin 
joints  are  located  arbitrarily.  The  bars  are  labeled  sequentially  from  0  to  3.  On  each  bar, 
a  body-fixed  frame  is  attached  such  that  its  origin  is  at  the  center  of  mass  of  the  bar  and 
the  x-axis  is  parallel  to  the  line  connecting  two  joints  on  the  bar.  In  particular,  we  choose 
the  positive  direction  of  of  z-axis  of  i-th  bar  towards  the  (i  +  1 ) -th  bar,  for  i  --  0, 1,2,3 
(mod  4).  We  define  the  following  quantities. 

d ij  the  vector  from  the  center  of  mass  of  i-th  bar  to  the  joint  with  j  -  th  bar 

in  body-fixed  frame; 

rj  the  position  vector  of  the  center  of  mass  of  i-th  bar  relative  to  an  inertial 

observer; 

r?  the  vector  from  the  system  center  of  mass  to  the  center  of  mass  of  i-th 

bar; 

rc  the  position  of  the  system  center  of  mass  relative  to  the  reference  point 

of  the  inertial  observer; 

Si  the  orientation  angle  of  i-th  bar  relative  to  the  inertial  frame; 

Sij  the  relative  angle  between  i-th  bar  and  j  -  th  bar,  i.e.,  SXJ  =  Si  —  Sj  ; 

li  the  length  of  i-th  bar,  which  is  defined  as  the  distance  between  the  joints 

on  i-th  bar,  i.e.  /,•  =  ||dt)j+1  -  d,)l-_1||; 

rrti ,  I{  the  mass  and  the  moment  of  inertia  of  i-th  bar  about  its  center  of  mass; 

m  the  total  mass  of  the  system,  i.e. 
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m  =  mi. 
i= o 

Any  pair  of  adjacent  bodies  is  connected  by  the  following  relation,  the  hinge 
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inertial 

observer 


Figure  2.1  The  general  structure  of  four-bar  linkage 


constraint , 


r?+i  =  r-  +  -  f?((9l+1)di+lil-  i  =  0,l,2,3  (mod  4).  (2.1) 


where 


R{9i)  = 


cos 


(&i)  —sin(di) 


sin(9i )  cos(9i ) 

the  rotation  matrix.  By  eliminating  r?  in  (2.1)  we  find  the  loop  constraint  or  closure 
constraint , 

3 

-d.,.-i)  =0,  (2.2) 


:=0 


where  the  convention  of  modulo  4  for  the  subscripts  is  adopted. 
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3.  THE  CONFIGURATION  SPACE 


In  this  section  we  investigate  the  conditions  under  which  the  loop  constraint  (2.2) 
describes  a  submanifold  of  the  configuration  manifold  of  an  open  four-bar  chain. 

For  a  planar  floating  four-bar  open  chain,  the  configuration  space  is 

M  =  R2  x  Sl  x  S1  x  S1  x  S1. 

Here  M  is  a  6-dimensional  smooth  manifold  with  local  coordinates 


q  =  (zo,  yo,  0oA,  #2,03)- 


This  corresponds  to  keeping  track  of  a  material  point  (say  center  of  mass  on  one  of  the 
bodies)  and  the  four  absolute  orientations.  See  (Oh  et  al.  [1989];  Sreenath  [1987];  Sreenath 
et  al.  [198S])  for  the  hamiltonian  mechanics  of  such  open  chains. 

For  a  closed  four-bar  mechanism  considered  in  this  paper,  the  configuration  space 
denoted  by  Q  is  a  subset  of  M  determined  by  (2.2),  i.e. 

Q  =  { q  €  M|F(q)  -  0} 


for 
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F(q)  =  £  K(0,  )(d.ii+1  -  d, 

1=0 


For  Q  to  be  a  submanifold  of  M ,  we  have  the  following  condition. 

Theorem  3.1:  If 


/o  i/l  ±b±l3  7^  0, 


(3.1) 


Q  is  a  submanifold  of  M . 

Proof.  Note  that  F  :  M  — >  R2 .  From  (Guillemin  &  Pollack  [1974])  we  know  that  if  0  is 
a  regular  value  of  F,  i.e.  OF/dq  has  full  rank  for  all  q  E  Q,  then  Q  is  a  submanifold  of 
M. 

From  (3.1)  we  have 

dF  __  f  0  0  —lo.sin(do)  —l\siri(9i)  —losin(9o)  —lisin(93)\ 

9q'\o  o  k  cos(9q )  licos{9i)  1 2 cos^do )  l^coslOi)  y 
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It  is  easy  to  check  that  all  the  nontrivial  determinants  of  2  x  2  submatrices  are  given  by 
the  following  functions, 

fib (q)  =  khsin(di  -  0O) 

02 (q)  =  lohsin(92  -  0O ) 

03 (q)  =  lohsin(93  -  0O)  (3.2) 

04(q)  =  hhsin(92  -  0a) 

05 (q)  =  l\hsin(d3  -  0i ) 

0e(q)  =  hhsin(93  -  d2). 

Therefore,  if  for  each  q  £  Q ,  there  exists  an  i  such  that  <7i(q)  ^  0  ,  Q  is  a  submanifold  of 
M .  It  is  obvious  that  the  above  condition  depends  on  the  relative  angles  and,  consequently, 
the  lengths  of  the  bars.  To  arrive  at  the  condition  in  the  statement  of  the  theorem,  we 
proceed  from  the  converse. 

If  0*(q)  =  0  for  all  i ,  from  (3.2),  we  have 

0i  —  0o  =  0  ov  7T  (3.3 a) 

and 

03  —  02  =  0  or  7r  (3.36) 

and 

$3  —  0O  =  0  or  7 r.  (3.3 c) 

Premultiplying  (2.2)  by  R(  —  9q),  we  get  following  equivalent  closure  constraint 
equations 

Zo  +  hcos(91  —  90 )  +  /2(cos(03  —  02)cos(03  —  0O )  +  sm(03  —  02)sm(03  —  0O)) 
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+  /,3  COS  (03  —  02  )  =  o 


(3.4a) 


Cases 

010 

032 

030 

Link  Conditions 

i 

0 

0 

0 

/0+  +  /2+  ^3=  0 

ii 

0 

0 

n 

lo+  h~  h~  0 

iii 

0 

7t 

0 

^o+  h~  h+  h=  0 

iv 

0 

K 

7C 

Iq+  ^1  ^2*  ^3=  ^ 

V 

7t 

0 

0 

V  ^1+  ^2+  ^3=  0 

vi 

K 

0 

TC 

/0-  /r  /2-  /3=  0 

vii 

K 

Jt 

0 

Iq-  /j-  /2+  /3=  0 

viii 

K 

n 

TC 

V  ^l+  V  ^3 =  0 

Structures 


a 


2_o L 


a 


-j-o - ^ - 

3  Q  0  o  1 


„Q  ,0  t  Q  2 


2  q  3  n  0 


1  Q  2  n  3 


3 _ n  0 

T~a - 1 — 

_ 3 


as 


0 


or 


0°— 1- 


XI 


s 


D 


D 


D 


D 

D 


Table  3.1 

l\sin(0i  —  #0)  +  h(cos(d 3  —  #2)-sin($3  —  #0)  ~  sin(6z  —  ^9)005(^3  —  #0)) 

+  Izsvn^Oz  —  $2)  =  0.  (3.46) 

So,  to  make  6{ ’s  satisfy  both  (3.3)  and  (3.4),  the  length  of  the  bars  should  satisfy 

6)  +  (  —  l)i'1h'r(— l)*”'2  L  +  (  — 1)^'3^3  —  0  (3-5) 


for  some  ^,^27^3  £  {0,1}.  By  contradiction,  if  for  all  ki,ko,k3  £  {0,1}  (3.5)  does  not 
hold,  then  there  is  an  i  such  that  </,-(q)  ^  0,  which  means  that  Q  is  a  submanifold  of  M . 


Remark  3.2:  Table  3.1  summarizes  the  conditions  on  lengths  of  the  bars  which  cause  Q 
not  to  be  a  submanifold  of  M  and  the  configurations  with  respect  to  the  relative  angles 
given  in  (3.3).  It  is  easy  to  observe  that  case  (i)  can  never  happen  since  lt  are  assumed  to 
be  positive.  In  addition,  cases  (iii),  (iv),  (v)  and  (vi)  are  trivial  since,  in  these  cases,  none 
of  the  relative  angles  can  vary,  i.e.  the  configuration  space  loses  one  degree  of  freedom  and 
the  linkage  becomes  a  rigid  structure.  | 

In  the  following  we  show  that  the  condition  in  Theorem  3.1  can  be  simplified  by 
ignoring  the  labels  on  the  bars  and  give  a  topological  description  of  Q.  We  first  recall 
some  definitions  and  results  in  the  classical  theory  of  mechanisms  (Grashof  [1883];  Paul 
[1979a];  Paul  [1979b]). 

For  a  four-bar  linkage  in  classical  mechanics,  in  which  one  bar,  the  frame  is  fixed, 
the  following  quantities  are  defined: 

5  =  length  of  shortest  bar 
/  =  length  of  longest  bar 
p,  q  —  lengths  of  intermediate  bars. 

A  bar  which  is  free  to  rotate  through  2?r  with  respect  to  a  second  bar  is  said  to  revolve 
relative  to  the  second  bar  and  is  referred  to  as  a  crank.  Any  bar  which  does  not  revolve  is 
called  a  rocker.  If  it  is  possible  for  all  bars  to  become  simultaneously  aligned,  such  a  state 
is  called  a  change  point  and  the  linkage  is  said  to  be  a  change-point  mechanism. 
Theorem  3.3:(Grashof  [1883]) 

(1)  A  four-bar  mechanism  has  at  least  one  revolving  link  if 

s  +  l  <  p  +  q 


and  all  three  will  rock  if 

s  +  l  >  p  +  q. 

(2)  A  four-bar  mechanism  is  a  change-point  mechanism  iff 


s  +  l  =  p  +  q. 
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Remark  3.4:  If  s  +  l  <  p  +  q,  the  shortest  bar  is  the  revolving  bar. 

It  is  easy  to  check  that  the  cases  (ii),  (vii)  and  (viii)  in  Table  3.1  correspond  to 
s  +  l  =  p  +  q,  i.e.  they  correspond  to  change-point  mechanisms.  From  Theorem  3.1,  we 
immediately  have  following  result. 

Corollary  3.5:  If  l  <  s  +  q  +  p  and 

s  +  l  ^p  +  q, 

Q  is  a  submanifold  of  M .  | 

Remarks  3.6: 

(a)  If  l  >  s  +  q  +  p,  the  mechanism  is  not  constructible. 

(b)  In  the  sake  of  convenience,  we  refer  to  as  the  condition  s  +  l  <  p  +  q  and  s  +  /  >  p  +  q 
as  Grashof  and  non-Grashof  condition,  respectively.  The  corresponding  linkages  are 
referred  to  as  Grashof  and  non-Grashof  mechanism,  respectively.  We  will  adopt 
the  same  terminology  for  floating  linkages  also.  So,  Corollary  3.5  says  that  Q  is  a 
submanifold  of  M  if  either  Grashof  or  non-Grashof  condition  holds. 

(c)  It  is  easy  to  see  that,  given  four  bars,  assembling  them  into  a  closed  loop  and  labeling 
them  sequentially,  the  linkage  may  have  instantaneously  one  of  three  kinds  of  shape. 
These  correspond  to 

f  >  0; 

sin(dz  -  d2)  <  <0; 

l  =  0. 

In  the  classical  theory  of  mechanisms,  they  are  called  lagging  form,  leading  form  and 
dead  point,  respectively.  g 

With  above  remarks  in  mind,  we  proceed  to  give  a  topological  description  of  the 
configuration  manifold  Q .  We  first  need  the  following  result. 

Proposition  3.7:  Let  l\  =  s.  Then  s  +  l  <  p  +  q  if  and  only  if  sm(03  —  62)  7-  0  for  all 
configurations. 

Proof.  The  mechanism  can  be  assembled  with  s  adjacent  to  /  or  with  s  opposite  l.  And, 
l  can  be  Iq,  I2  or  I3  .  If  O3  —  62  —  kiv ,  the  whole  structure  attains  a  triangular  shape  which 
has  the  property  that  the  sum  of  two  sides  is  larger  than  the  third  one.  Then  it  is  easy  to 
check  that  all  possible  cases  will  lead  to 

s  +  l  >  p  +  q. 
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Further,  it  is  obvious  that  if  s  +  1  =  p  +  q,  there  exist  #2  and  83  such  that  sm(03  —  02)  =  0. 
And,  if  5  +  l  >  P  +  q,  by  Grashof’s  theorem,  all  bars  will  rock  with  respect  to  each 
other.  This  means  that  dead  point  is  reachable.  Therefore,  if  sin(8 3  —62)  f  0  for  all 
configurations,  s  +  l  <  p  +  q.  g 

Remark  3.8:  An  equivalent  way  to  state  above  assertion  is  that  a  four-bar  linkage  is  a 
Grashof  mechanism  if  and  only  if  it  is  constrained  to  be  in  either  leading  form  or  lagging 
form.  Moreover,  for  non-Grashof  mechanisms,  the  linkage  can  vary  continuously  from 
leading  form  to  lagging  form.  g 

From  Grashof’s  theorem  and  the  above  proposition  we  can  get  a  topological  descrip¬ 
tion  for  Q . 

Theorem  3.9: 

(a)  For  a  Grashof  linkage,  i.e.  s  +  l  <  p  +  q,  Q  has  two  components.  Each  component 
is  diffeomorphic  to  R2  x  S1  X  S1 . 

(b)  For  a  non-Grashof  linkage,  i.e.  s  +  l  >  p  +  q,  Q  —  R2  x  51  x  S1. 

Proof.  Our  proof  is  based  on  explicit  parameterization  of  the  manifold  Q.  Recall  that 
from  (3.1)  the  dimension  of  Q  is  four.  Again,  we  let  s  =  l\ . 

(a)  If  s  +  /  <  p  +  q,  we  consider  the  parameters  (xo,  yo,  #0, 8\ ) ,  where  ( Xo,yo )  is  the 
coordinate  of  any  point  on  0-th  bar  in  inertial  frame.  From  the  definitions  of  80  and  8y 
and  the  Grashof  theorem,  both  8q  and  8 1  can  vary  from  —  7r  to  tt  independently.  From 
the  remark  following  Proposition  3.7,  one  component  of  the  manifold  Q  corresponds  to 
leading  form.  The  other  corresponds  to  lagging  form. 

(b)  If  s  +  l  >  p  +  q,  from  Grashof  theorem,  there  exists  an  angle  a ,  0  <  a  <  n ,  such 
that,  #10  G  [— a  +  7r,  a  +  7r] .  At  the  boundaries,  the  system  is  at  “dead  point”  (see  Remark 
3.6  (c)).  Now,  we  consider  the  independent  parameters  (xo,  yo ,  8q,  6) ,  where 

f ,  if  sm(032)  >  0  (lagging  form); 

,  if  sin(9 32)  <  0  (leading  form). 

It  is  easy  to  check  that  9  and  #10  have  a  one-to-one  relation  in  both  leading  form  and 
lagging  form.  Moreover,  9  varies  from  —  tv  to  7 r.  For  8  G  (— 7r,0),  the  mechanism  is 
in  the  leading  form;  for  9  G  [0,7r],  the  mechanism  is  in  lagging  form.  So  (xq  ,  yo ,  $0 , 9) 
parameterize  R2  X  S1  X  S1  .  g 
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Remarks  3.10: 

(1)  Gibson  and  Newstead  in  (Gibson  &  Newstead  [1986])  proved  above  assertion  through 
the  methods  of  algebraic  geometry. 

(2)  This  theorem  says  that  for  fixed  point  in  R 2  ,  the  configuration  space  is  the  disjoint 
union  of  two  tori  for  a  Grashof  linkage  and  a  torus  for  a  non-Grashof  linkage.  For 
the  latter,  the  torus  can  be  split  into  two  parts,  one  corresponding  to  leading  form, 
the  other  corresponding  to  lagging  form. 

(3)  From  the  above  remark,  we  note  that,  for  both  Grashof  and  non-Grashof  linkages, 

one  may  use  (xo,  j/o,  9q,  )  and  specified  form  (leading  or  lagging)  to  parameterize 

Q  locally.  For  non-Grashof  case,  one  has  to  worry  about  the  parameterization  in  the 
neighborhood  of  dead  point.  This  problem  can  be  solved  by  re-labeling  the  bars.  We 
will  discuss  it  further  in  the  next  section.  | 


Ox 


Figure  3.1  The  reduced  configurations,  an  example. 
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Example  3.11:  Here  we  illustrate  the  Theorems  3.1  and  3.9  by  a  simple  example.  First, 
we  fix  one  bar  (position  and  orientation),  say  Iq  and  let  #0  =  0.  This  means  a  point  in 
R2  x  Sl  has  been  chosen  and  the  dimension  of  the  configuration  is  reduced  to  one.  Now 
the  constraint  equation  (2.2)  becomes 

lo  +  IxCOS^d i)  T  /2co5(^2)  T  ^3 cos(^3 )  =  0 

lisin(di)  +  l2sin(02)  +  hsin(6z)  —  0. 

Eliminating  one  more  angle,  say  9 3  ,  we  get 

/($i ,  $2)  =  (lo  T  l\cos(^9\ )  +  /2co‘S(^2 ))2  +  (l\sin(9i )  4"  /2sin.($2 ))2  —  l\  —  0.  (2-6) 

The  solutions  of  this  equation  gives  a  curve  on  a  torus  T2  as  the  configuration  space. 

Now  we  choose  Iq  =  3,  l\  =  3,  I3  =  4  and  let  I2  vary.  Figure  3.1  shows  the 

results.  The  rectangle  with  opposite  edges  identified  is  the  standard  way  to  represent  a 

torus  (Guillemin  &  Pollack  [1974]).  We  see  that  when  I2  —  1.9,  s  +  l  <  p  +  q  holds  and 
the  solution  of  (3.6)  on  T 2  is  a  disconnected  closed  curve;  when  I2  =  2,  s  +  /  =  p  +  q 
holds  and  the  solution  of  (3.6)  on  T2  is  a  “figure  8”;  when  I2  =  3,  s  +  l  >  p  +  q  holds  and 
the  solution  of  (3.6)  on  T2  is  a  connected  closed  curve;  when  /2  =  4,  s  +  l  =  p  +  q  holds 
and  the  solution  of  (3.6)  on  T 1  is  a  “figure  8”,  again;  when  /2  =  4.5,  s  +  l  >  p  +  q  holds 

and  the  solution  of  (3.6)  on  T2  is  a  single  closed  curve,  again.  If  the  linkage  is  allowed  to 

float,  we  get  the  configuration  spaces  described  in  Theorem  3.9.  | 
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4.  KINETIC  ENERGY 


In  this  section  we  derive  the  kinetic  energy,  or  lagrangian  since  we  assumed  that  no 
potential  energy  is  involved,  for  the  whole  system.  The  basic  idea  is  to  write  the  kinetic 
energy  for  each  individual  body  first  and  then  use  the  constraint  equations  to  eliminate 
extra  variables. 

The  kinetic  energy  of  the  i-th  bar  is 


1  .2  T.  ,  1  . Ili.l|2 


Ti  =  -ut  It  +  -mi||r,- 
where  u>i  =  9t .  The  total  kinetic  energy  is 

3  3 

T  =  9  ^  0  X/  miH^I|2- 

1=0  2=0 

To  describe  the  kinetic  energy  relative  to  the  center  of  mass,  we  have  following  useful 
equations, 


r;  =  rc  +  r?  *  =  0, 1,2,3 

3 

rnlrci  =  0. 
i=0 

By  applying  (4.2)  and  (4.3),  (4.1)  becomes, 

3  3 

T  =  +  +  ^INI2- 

1=0  2=0 

Applying  (2.1)  and  (4.2),  we  get 

r  i  =  — 
m 

-  RidiXmi-idi'i-i  +  {mi+i  +  ml+2)di,i+i) 

+  R(^1+1  )((mi+ 1  +  mi+2)di+i)i  —  mi+2di+i)i+2) 


(4.2) 

(4.3) 


(4.4) 


+  R(6,i+2)mi+2d!+2Ii+i] 

for  i  =  0,1,2, 3  (mod  4).  Note  that  the  convention  on  the  subscript  can  be  used  because 
of  the  closed-loop  condition.  Furthermore, 

r  •  =  —  {mi-iLbi-iR(dl-i)di-itl 

rn 

-  d)iR(0i)(7^i-idi,i'-i  +  (mi+ 1  +  mi+2)dl,J+i) 

+  Wj+ii?(^;.|.i)((m;+i  4-  rn.;_(_2)d;+iji  —  m;4_2di4-i  ,i+2) 


+  l^i+2R(9i+2)mi+2di+2,i+l] 


(4.5) 
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for  i  —  0, 1, 2,  3  (mod  4),  where 


By  substituting  the  formula  for  r ?  into  (4.4),  we  get  a  more  compact  expression  for  the 
kinetic  energy 


T  =  ^  <  u,  Ju  >  +im||rc||2  (4.6) 

where  u  —  (tao ,  o>x ,  u>2 ,  <^3  )T  and  J  =  ( =  0,1,  2, 3)  is  a  4  x  4  symmetric  matrix, 
with  elements  given  as  follows. 

Let 


M-  -  —  [mtm,+i(m,-  +  m,+i) 
m *■ 

+  mi+i(m1+1mi_i  +  m;mi+ 2) 

+  mtmi+2(mi+1  +  mt+2 )]  (4.7a) 

Mj1  =  ~  (4.76) 

(4.7c) 

Then 


Jit  =  T  +  M/lldll^^M^lldll^, 

—  2 M/i1  <  d^j+i,  di^-i  >  (4.8a) 

^1,1+1  ~  d j ; 1 , 77(^i-)-i , i ) d i ^  > 

+  Ml1  <  d  t.i+l,  -R(^i+l,i)di+l,i+2  > 

+  M/_f1  <  d,'^-!,  i?(0i+1)i)d,+  lji  > 

+  M-11  <  d^i-j ,  i?(0i+1),)d,+li;+2  >  (4.86) 

J i.t+2  —  —Mj1  <  d,  i+i,  i?(0i+2,i)di+2>i+l  > 

—  M/_^(  <  d^j+x, i?(5l+2,i)dj+2,i-i  > 

—  Mj_l_ 2  <  d^x  —  i ,  i?(0;+2,i)d;_|_2,i  —  1  > 

—  Mf11  <  d^i-i,  i?.(0i+2,i)dt+2,i+i  >  (4.8c) 
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for  i  —  0, 1, 2, 3  (mod  4). 

Remark  4.1:  In  (4.6),  8i  and  the  velocities  uq,  i  =  0,1, 2, 3  are  involved.  Since  the 
whole  linkage  forms  a  closed  loop,  these  variables  are  not  independent.  Let  l\  —  s . 
By  the  Implicit  Function  Theorem  one  can  show  that  given  a  point  (80 ,  ,  02,  83)  in 

S1  x  51  x  S1  x  51 ,  such  that  (2.2)  and  sin(6z  —  62)  7^  0  are  satisfied,  around  this  point 
there  exists  a  neighborhood  in  which  81  and  8 3  can  be  represented  as  the  functions  of 
(8q,  8\)  uniquely.  For  a  Grashof  linkage,  from  Proposition  3.7,  this  condition  is  satisfied 
everywhere  in  Q .  But  for  non- Grashof  linkages  this  condition  is  satisfied  almost  everywhere 
in  Q ,  except  at  dead  points.  We  note  however  that  even  in  this  case  all  four  bars  cannot 
be  aligned.  So  at  such  points,  8q  and  6\  can  be  represented  as  the  functions  of  (82,83) 
uniquely.  This  is  equivalent  to  relabeling  the  four  bars  such  that  sin(83  —  82)  ^  0. 
Moreover,  from  (Paul  [1979a] ),  locally,  there  exist  two  functions  f\  and  /2  of  (8 1  —  8q) 
such  that 

82  —  8q  +  /1  (8\  —  8q  ) 


83  —  8q  +  /2($i  —  0o). 


We  note  that  the  expressions  for  functions  f\  and  /2  depend  on  the  sign  of  sin(^3  —80), 
i.e.  on  whether  the  linkage  is  in  leading  form  or  lagging  form.  Therefore,  we  conclude 
that  at  every  point  on  Q ,  every  element  of  the  matrix  J  can  be  expressed  as  a  function 
of  8\  $0 .  | 


Differentiating  the  loop  constraints  (2.2),  one  can  get  a  relation  between  (cuo,u;i) 


and  (w2,w3): 


(4.9a) 


where 


(4.96) 

\  tasin(os-t>2)  i3sm(o3-o2)  / 

Again,  following  the  same  line  of  reasoning  as  in  Remark  4.1,  the  matrix  D  is  well-defined 
locally  in  general. 

We  summarize  the  above  discussion  in  the  following  theorem. 


n  = 


__  J  I  V3  ; 

l2sin(63  —  d2) 
l0sin(82-60) 


/2S!7l(53  —  (?2  ) 
l,sin(92  —  d  1 ) 
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Theorem  4.2:  The  kinetic  energy  of  a  floating  four-bar  linkage  can  be  represented  as 


T  =  i  <  <jj,  Jw  >  +^m||rc||2  (4.10) 

where  ui  =  and 

j  =  d  «nj(') 

for  J  given  in  (4.6)-(4.8)  and  fl  given  in  (4.9).  The  elements  of  J  are,  locally,  branched 
functions  of  (9\  —  d0)  (depending  sgn(sin(0s  —  62))).  I 

Before  ending  this  section,  we  give  a  property  of  the  matrix  which  will  be  used  in 
Sections  6  and  7. 

Proposition  4.3:  If  fl  is  well  defined, 


Proof:  Premultiplying  (2.2)  by  R(— 92)  and  R(—9, 3),  we  get 

losin(0o  —  62)  +  l\sin(9i  —  92)  +  hsin{9j,  —  92)  =  0 


and 

losin(90  —  #3)  +  l-[sin(9i  —  $3)  +  l2sin(02  —  #3)  =  0 
respectively.  The  result  follows  immediately. 


■ 
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5.  SYMMETRY,  INTEGRAL  AND  REDUCTION 


We  shall  show  here  that  a  floating  four-bar  linkage  is  a  simple  mechanical  system 
with  symmetry  in  the  sense  of  Smale  (Abraham  &  Marsden  [1978];  Smale  [1970]). 

A  simple  mechanical  system  with  symmetry  is  a  4-tuple  (Q,  K,  V,  G),  where, 

(i)  (Q,K)  is  a  Riemannian  configuration  manifold  with  metric  K ; 

(ii)  G  is  a  Lie  group  acting  on  Q  on  the  left, 

<f>  :  G  x  Q  —*  Q 

(</,?)  >-»  $<»(?)  =  $(9,<l) 

such  that  for  each  g  €  (7,  is  an  isometry  of  (Q,  Jv); 

(iii)  V  :  Q  — >  R  is  a  G-invariant  potential  function.  The  associated  lagrangian  is 
defined  by 


1 


(5.1) 


L  :  TQ  ->  R 

vq  ^  L{vq)  =  ^ K(vq,vq )  -Vo  r(vg) 

where  r  :  TQ  —>  Q  is  the  canonical  tangent  projection.  The  Legendre  transform  FL  of  L 
is  given  here  by  the  vector  bundle  isomorphism 


IC  :  TQ  -+  T*Q 


satisfying 


Ka{Vg)  ■  Wg  =  K(Vg,Wg)  V  V  g  ,  W  g  6  TgQ. 


The  hamiltonian  H  :  T*Q  — >  R  is  defined  to  be 


H(a,)  =  i/i:((A'‘)-1(a,),(A-k)-1(a,))  +  r<.T'(a,)  (5.2) 

where  t*  :  T*Q  —>  Q  is  the  canonical  cotangent  projection.  Letting  Ro  denote  the 
canonical  symplectic  structure  (Abraham  &  Marsden  [1978])  on  T*Q,  the  hamiltonian 
dynamics  on  T*Q  is  given  by  the  unique  vector  field  Xh  on  T*Q  such  that 

dH(Y)  =  n0(XH,Y) 

for  any  vector  field  Y  on  T*Q . 

IS 


Let  :  G  X  T*Q  — >  T*Q  be  the  cotangent  lift  of  the  action  <f>.  Denoting 

*ro = we  have 

H  o  =  JET,  .  (5.3) 

i.e.  the  group  G  is  a  symmetry  group  of  the  hamiltonian  system  (T*Q,Q0^Xh) • 

The  modern  setting  of  Nother’s  theorem  relating  symmetry  to  the  existence  of 
integrals  of  motion  is  given  by  the  concept  of  momentum  mapping.  Let  Q  denote  the 
Lie  algebra  of  G  and  Q*  its  dual.  The  map 

J  :T*Q^  Q* 

given  by 

J(a,)  •(=«,(,(()  Vo,Ei;g  (5.4) 

is  Ad* -equi variant,  where  £  G  Q  and  £q  is  the  infinitesimal  generator  of  <f>  on  Q  associated 
to  £  (see  (Abraham  &;  Marsden  [1978])  corollary  4.2.11).  J  is  a  momentum  mapping  and 
the  G-invariance  of  the  hamiltonian  H  implies  that  J  is  an  integral,  i.e.  it  is  conserved 
along  trajectories  of  Xh  • 

Returning  to  four-bar  linkages,  we  restrict  attention  for  the  moment  to  linkages  of 
Grashof  type,  i.e.  the  condition  s  +  l  <  p  +  q  holds.  For  non-Grashof  case  we  will  give  a 
discussion  in  the  final  remark  of  this  section.  To  fix  the  parameterization,  we  let  l\  =  s .  We 
also  restrict  ourselves  to  a  connected  component  of  the  configuration  space,  i.e.  we  assume 
the  linkage  has  been  in  either  leading  or  lagging  form.  As  we  have  seen  in  previous  sections, 
the  lagrangian  only  depends  on  90  —  9\  and  9q  and  9\  .  Therefore,  in  some  essential  ways, 
the  description  of  the  system  will  parallel  that  of  a  coupled  two-body  problem  (Sreenath 
[1987]).  In  the  following,  we  will  skip  some  details  of  the  proofs. 

Note  that  the  lagrangian  in  (4.10)  does  not  depend  on  the  location  of  the  center  of 
mass  of  the  system.  This  means  that  the  dynamics  of  the  system  is  invariant  under  the 
translation  in  inertial  space.  Without  loss  of  generality,  one  can  place  the  inertial  observer 
at  the  center  of  mass  of  the  system.  In  (Sreenath  et  al.  [1988])  this  process  is  explained  via 
symplectic  reduction  by  the  translation  group  R 2 .  With  this  choice  of  inertial  observer, 
the  configuration  space  becomes 

Q  =  S1  X  S1  (5.5) 
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and  the  kinetic  energy  becomes 

T  =  -  <  uj,  Ju>  >  . 

The  corresponding  kinetic  energy  metric  on  Q  can  be  defined  as 

K(w i,w2)=  -  <wi,Jw2>  (5-6) 

for  iwi,  u>2  E  TqQ .  It  is  easy  to  observe  that  the  group  of  rigid  motions  of  this  system 
is  the  symmetry  group  G  =  Sl .  Denoting  a  point  in  Q  by  (0o,  9 j),  the  action  $  of  this 
group  on  Q  is  given  by 

(@o,  9\ ))  =  ($o  +  <j>,  +  <!>)•  (5-7) 

The  hamiltonian  is  given  by 

H  =  \<P,  J“V>  (5-8) 

for  n  =  Kh(u>)  =  (no,fii).  Using  the  abstract  formula  (5.4)  and  the  action  (5.7),  one  can 
show  that  a  momentum  mapping  for  the  action  $  is  represented  by 

•  £  =  no  +  Hi  (5.9) 

where 

«-(! 

Of  course,  u  =  /i0  +  Hi  is  conserved  along  trajectories  of  Xh  for  the  hamiltonian  H  in 
(5.8)  and  it  is  simply  the  net  angular  momentum  of  the  floating  four  bar  linkage  relative 
to  an  observer  at  the  system  center  of  mass. 

The  dynamical  trajectories  are  confined  to  a  level  set  of  the  form  J-1(z/).  The 
group  Sl  viewed  as  the  isotropy  subgroup  of  the  momentum  value  acts  freely  on 

and  one  gets  the  symplectically  reduces  dynamics  Xh„  on  the  reduced  phase 
space  Pu  =  J~1(u)/S1  ~  S1  x  R1  . 

As  in  (Sreenath  et  al.  [1988])  it  is  also  possible  to  Poisson-reduce  the  dynamics.  We 
recall  that  given  a  symplectic  manifold  (M,  u>),  and  a  smooth,  free,  proper,  symplectic 
action  of  Lie  group  G  on  M ,  the  canonical  Poisson  structure  on  M  defined  by 

{f,g}M=u{Xf,Xg)  V/,  g  e  C°°(M) 
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descends  to  a  Poisson  structure  on  the  quotient  P  =  M/G .  The  latter  is  defined  by 

{AsJm/gO’t  =  {/o  tt,<7  o  tt}m  (5.10) 

where/, g  £  C°°(M/G )  and  tt  :  M  — »  M/G  is  the  canonical  projection.  If  H  :  M  — >  i?  is 
a  G-invariant  hamiltonian,  it  induces  H  :  M/G  R  defined  by  H  o  tc  =  H .  We  refer  to 
the  dynamics  (vector  field)  X ^  defined  by 

Xf,(f)  =  {f,H}  VfeC™(M/G)  (5.11) 


as  the  Poisson  reduction  of  the  dynamics  X h  • 

In  the  present  context,  with  Q,K,H  as  in  (5.5)-(5.8),  the  space  M  —  T*(SX  x  S1) 
with  parameterization  (0q,  Q\ ,  po,  Mi)  carries  the  Poisson  structure, 


{/,*}  =  Bf 


!  =  0 


dg 

dm 


dl'dg_ 

dm  ddi 


) 


(5.12) 


for  all  /,  g  €  C°°(T*(5'1  x  S1)).  The  action  of  G  —  Sl  on  Q  given  by  (5.7)  is  free  and 
proper.  The  quotient  P  =  T*(Sl  x  S1)/Sl  ~  51  x  R2  carries  a  reduced  Poisson  structure. 
Parameterizing  P  =  T*(5'  x  Sl)/Sl  by  (01O,  m  ■>  m  )  >  the  Poisson  bracket  on  P  is  given 


by, 


{/>£} 


T* (Sl  xSl)/Sl  = 


df  ,  d<7 


(• 


dg 


df  df 


dd  io  5/ii  <9po  <9$io  d/i  i  <9po 

which  is  a  noncanonical  structure.  The  reduced  hamiltonian  H  is  given  by 


) 


(5.13) 


H (0io,  ^o,  Hi)  —  m  )^ 


(5.14) 


since  J  in  (5.8)  is  a  matrix  of  functions  of  the  difference  9\q  =  9\  —  9q  only.  The  reduced 
dynamics  is  then  immediately  given: 


Ao 


dH 

d0  io 


m  = 


dH 
dd  io 


dH 

dm 


dH 

df-io 


(5.15) 
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Equation  (5.15)  involves  complicated  analytic  expressions  resulting  from  the  substitutions 
for  #3  and  #2  in  terms  of  90  and  6 i  as  in  (Paul  [1979a]).  Certain  qualitative  aspects  of  the 
reduced  dynamics  can  still  be  explored  such  as  relative  equilibria  which  we  will  investigate 
in  the  following  section. 

Remark  5.1:  For  non-Grashof  mechanisms,  the  system  still  has  symmetry  group  G  —  Sl . 
This  is  clear  if  one  thinks  a  four-bar  linkage  as  an  open  four-bar  chain  combined  with 
constraint  equation  which  is  invariant  under  the  action  of  G .  It  is  also  clear  that,  in  this 
case,  if  #o  and  6\  are  just  parameters  of  Sl  X  S1,  instead  of  real  angles  of  0-th  and  1-st 
bars,  all  of  results  in  this  section  are  valid.  However,  for  convenience,  in  next  two  sections 
we  will  still  use  physical  angles  dQ  and  6 i  for  non-Grashof  mechanisms.  Then,#,)  and  6\ 
will  play  the  role  of  local  coordinates  for  S 1  x  S 1  and,  consequently,  (5.9),  (5.13)  and  (5.15) 
will  be  the  expressions  of  momentum,  reduced  Poisson  bracket  and  reduced  dynamics  in 
local  coordinates,  respectively.  | 

One  of  the  advantages  of  applying  reduction  theory  to  mechanical  systems  is  that 
it  helps  to  make  the  dynamics  of  the  system  more  transparent.  For  our  system  which  is 
of  four  dimensions  in  phase  space,  the  reduction  process  make  it  possible  to  display  the 
dynamics  by  the  phase  portrait  on  a  two-dimensional  surface.  To  illustrate  this,  we  show 
two  examples  here. 

Example  5.2: 

Consider  a  floating  four-bar  linkage  whose  parameters  satisfy  the  Grashof  condition 
and  which  is  of  lagging  form.  In  particular,  the  parameters  are  chosen  as  the  following. 

m0  =  1,  mi  =  1,  m2  =  1,  m3  =  1; 

i"o  ==  1 1  /i  —  1 ,  /2  —  1  j  I's  —  1 5 

d03  =  (-1.5,1),  doi  =(1.5,1),  dlo  =  (-0.5,1.3),  d12  =  (0.5, 1.3), 

d21  =(-1.5,1),  d23  =  (1.5,1),  d32  =  (—2,  A),  d30  =  (2,  A). 

Following  the  same  procedure  as  in  (Sreenath  et  al.  [19SS]),  the  dynamics  can  be 
further  reduced  to  a  symplectic  leaf.  The  hamiltonian  on  a  leaf  is  a  function  of  #10  and 
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v(=  \{}H  —  Ho))  and  the  dynamics  on  the  leaf  is  given  by 


d9 10 
dt 


9H 

lb 


(9io,v), 


dv 

dt 


9H 
99 10 


(010, 


Therefore,  for  fixed  value  of  angular  momentum,  one  can  draw  the  phase  portrait  on  a 
cylinder  with  #io  versus  u.  Figure  5.1  shows  the  phase  portraits  for  A  =  0  and  A  =  —5, 
respectively.  On  Figure  5.1(a),  we  see  there  is  one  center  and  one  saddle  point.  However, 
in  Figure  5.1(b),  there  are  two  centers  and  two  saddles.  The  difference  is  caused  by  an 
offset  in  the  position  of  center  of  mass  of  3rd  bar  from  the  line  connecting  the  two  joints 
on  that  bar.  | 

Example  5.3: 

Consider  a  floating  four-bar  linkage  whose  parameters  satisfy  the  non-Grashof 
condition  and  the  parameters  are  chosen  as  the  following. 


mo  =  1,  m  i 

=  1, 

m2  =  1,  m3  =  1; 

h  =  1,  I\ 

=  1, 

r— 1 

II 

i — f 

II 

do3  =  (-2,  A), 

CL 

o 

II 

JO 

1 

lcT 

l — ♦ 

II 

o 

d12  =  (1.5,-1), 

d21  =  (-1.5, -1.4),  d23  =(1.5, -1.4),  d32  =  (-1.5, -1),  d30  =  (1.5, -1). 

Instead  of  displaying  the  dynamics  on  position-momentum  phase  space,  here  we  try  to 
show  it  on  position- velocity  phase  space.  Since  the  dynamics  cannot  be  further  reduced  so 
that  it  depends  on  one  relative  angle  and  corresponding  angular  velocity,  one  has  to  display 
the  phase  portrait  in  three  dimensional  space.  For  a  given  angular  momentum,  the  phase 
portrait  sits  on  a  surface  in  this  space.  In  this  example,  this  space  is  parameterized  by 
(dioHos)  -  Figure  5.2  shows  the  phase  portraits  corresponding  to  A  =  —15  and  A  =  0, 
respectively.  From  Figure  5.2,  we  see,  again,  the  change  of  numbers  of  centers  and  saddles. 

I 

In  the  next  section,  we  will  show  how  to  compute  the  centers  and  saddle  points  and 
associated  bifurcations. 
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6,  RELATIVE  EQUILIBRIA  AND  BIFURCATIONS 

We  first  give  the  definition  of  relative  equilibria  and  recall  Smale’s  theorem  (Smale 

[1970]). 

For  a  hamiltonian  system  (M, uj,Xh)  where  M  is  a  2n -dimensional  smooth  mani¬ 
fold,  u  is  a  symplectic  form  on  M ,  Xh  is  hamiltonian  vector  field  determined  by 

ixHu  =  dH 

for  some  smooth  hamiltonian  function  H  on  M.  For  simple  mechanical  systems, 
M  =  T*Q,  the  cotangent  bundle  of  the  configuration  space  Q,  of  the  system.  Let  G 
be  a  symmetry  group  acting  on  Q  on  the  left, 

$  :  G  x  Q  -+  Q 

(g,q)  ^  $s(?)  =  Hg,q), 

and  hence  on  T*Q  by  cotangent  lift.  The  following  definition  of  a  relative  equilibrium  is 
standard  (Abraham  &  Marsden  [1978]). 

Definition  6.1:  Let  Fy  be  the  flow  of  Xh  on  M .  Then  2e  G  M  is  a  relative  equilibrium 
if  Fxu(ze)  is  a  stationary  motion,  i.e.  there  exists  £  £  Q  such  that 

F!\aM  =  eXP(f  0(-e), 

where  Q  is  the  Lie  algebra  of  the  group  G .  | 

Remarks  6.2: 

(1)  Let  Xfj  be  the  Poisson  reduced  vector  field  as  shown  in  (5.11).  Then  3e  is  a  relative 
equilibrium  iff 

4W-7))  =  o. 

(2)  A  physical  interpretation  of  relative  equilibria  is  that  if  the  dynamics  of  a  system  is 
rotationally  invariant,  the  dynamical  orbit  of  a  relative  equilibrium  appears  to  be  a 
fixed  point  for  an  observer  in  a  suitable  uniformly  rotating  coordinate  system.  | 
Given  a  simple  mechanical  system  with  symmetry,  (Q,  I\,  V,  (7),  which  is  defined  in 

section  5,  we  have  following  theorem  by  Smale. 
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Hence,  by  Theorem  6.3,  (0e,/ie)  is  a  relative  equilibrium  point  on  T*Q  iff  9e  is  a  critical 
point  of  the  function 

where  9e  =  ((#o)e,  ($i)e)  and  =  ((po)e,  (pi)e)-  Applying  Proposition  4.3  and  the  4  x  4 
matrix  J  defined  in  Theorem  4.2,  we  have 

^(<Mi)  =  -£2eTJe  (6.5) 

where  e  =  (1  1  1  1)T.  Since  the  elements  of  matrix  J  are,  locally,  functions  of  9\  —  90  or 
constants,  the  above  V£  satisfies  (6.2)  for  all  g  6  Sl .  It  follows  that 

y?(M  =  n(0o,0i). 

Then,  the  critical  points  9e  of  will  make  ((#i)e  —  (#o)e)  to  be  the  critical  points  of  . 
At  relative  equilibrium,  the  relative  angles  between  bars  keep  unchanged  and  the  whole 
system  rotates  around  the  system  center  of  mass  with  constant  angular  velocity. 

Remarks  6.5: 

(1)  It  should  be  noted  that,  in  general,  given  the  value  of  ((#i)e  —  (#o)e)  one  cannot  tell 
what  the  relative  equilibrium  shape  of  the  linkage  looks  like.  The  particular  form 
(leading  or  lagging)  has  to  be  indicated  at  the  same  time. 

(2)  with  £  =  1  is  the  locked  inertia.  Its  value  at  @10  equals  to  the  value  of  the  moment 

of  inertia  of  the  corresponding  frozen  system,  i.e.  the  system  with  all  joints  locked, 
about  its  center  of  mass.  The  above  result  shows  that  a  relative  equilibrium  shape 
corresponds  to  a  frozen  system  which  has  maximum  or  minimum  value  of  moment 
of  inertia  within  all  possible  frozen  systems.  Since  at  relative  equilibrium  if  is  a 
constant,  we  will  let  £  =  1  later.  | 

In  the  rest  of  this  section,  we  are  particularly  interested  in  assemblies  which  admit 

configurations  with  reflection  symmetry,  which  will  be  called  symmetric  configuration. 
Applying  the  notations  in  Section  2,  a  floating  four-bar  linkage  is  of  symmetric  type  if, 
with  proper  consecutive  labeling  of  the  bars, 

mi  =  m3,  lx  =  /3  (6.6a) 
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(6.6 b) 


|doi|  =  |do3 1,  |dio|  =  |d30|,  jdi2 1  =  |d32|,  ]d2i|  =  jd23 

and 

df0  (q  d30  >  0.  (6.6c) 

In  other  words  it  can  form  symmetric  shapes  as  shown  in  Figure  6.1.  It  is  not  hard  to  verify 
that  any  four-bar  linkage  which  satisfies  (6.6)  has  two  such  symmetric  configurations. 


Figure  6.1  A  Symmetric  Configuration 

Although  it  is  not  easy  to  find  the  critical  points  of  analytically,  for  a  particular 
example  one  can  easily  find  them  numerically.  Unlike  the  planar  two  body  problem 
(Sreenath  [1987])  for  which  the  dimension  of  the  shape  space  is  also  one,  here  the  function 
Vj  has  many  parameters.  A  natural  question  is  to  determine  how  these  parameters 
affect  the  relative  equilibria,  e.g.  their  numbers  and  location  on  shape  space,  etc..  Of 
course,  it  is  difficult  to  answer  this  question  for  completely  arbitrary  choice  of  parameters. 
However,  by  leaving  one  particular  parameter  free  such  that  the  assembly  preserves 
its  symmetric  configurations  aird  freezing  all  other  parameters,  one  still  can  observe  a 
nontrivial  bifurcation  phenomenon.  To  illustrate  this  we  consider  an  example. 
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Example  6.6:  Let  us  choose  the  parameters  as  follows. 


rn0  =  1,  mj  =  1,  m2  =  1,  m3  =  1; 

do3  =  (— 2,Ao),  doi=(2,  Ao),  dio  =  (— 1.5,  — 1),  di2  =  (1.5,  — 1), 

d21  =  (-1.5,  -1.4),  d23  =  (1.5,  -1.4),  d32  =  (-1.5,  -1),  d30  =  (1.5,  -1). 

Now  the  assembly  has  non- Grashof  structure. 

Using  Uj ,  for  any  Ao  one  can  find  relative  equilibria  (#io)e  for  both  leading  form 
and  lagging  form,  and  hence,  corresponding  9  which  is  defined  in  the  proof  the  Theorem 
3.9.  As  Ao  varies  from  — oo  to  +oo,  one  can  plot  a  diagram  for  9  versus  Ao .  Figure 
6.2  shows  the  result,  in  which  solid  dots  represent  stable  relative  equilibria,  small  circles 
represent  unstable  relative  equilibria.  | 


Figure  6.2  Bifurcation  diagram:  an  example 
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From  this  example  one  can  make  the  following  empirical  observations: 

(1)  There  are  two  unbounded  symmetric  branches  on  the  diagram  and  and  these 
branches  are  bifurcated  at  some  points.  The  bifurcations  appear  to  be  pitchfork  bifurca¬ 
tions. 

(2)  Almost  any  value  of  9  can  be  a  relative  equilibrium  for  a  particular  Ao  .  In  other 
words  the  bifurcation  diagram  is  connected  globally. 

(3)  The  number  of  relative  equilibria  can  be  two,  six  and  ten. 

The  first  observation,  which  relates  to  the  local  bifurcation  problem,  is  what  we  will 
concentrate  on  in  the  rest  of  this  section.  The  others  will  be  discussed  later. 

As  we  have  seen,  the  function  of  a  four-bar  linkage  is  a  multiple  parameter 
function.  One  might  expect  very  complicated  bifurcation  features  with  respect  to  these 
parameters.  Here,  instead  of  considering  a  general  structure,  we  study  a  special  assembly 
which  has  symmetric  configuration.  To  avoid  too  many  tedious  calculations  we  particularly 
choose  the  parameters  of  the  assembly  as  follows. 


m0  —  mi  =  m2  =  rn$  —  1; 


=  (—do,  A0), 

doi  —  (do,  A0), 

dio  —  (—TO), 

di2  =  (l,0) 

(6.7) 

=  (-d2,  A2), 

d23  =  ( do ,  A2), 

d32  =  (  —  1, 0), 

d.30  =  (1,0). 

where  do  and  d2  are  fixed  and  do  >  do  >  0,  Ao,A2  G  R .  Moreover,  we  consider  the 
non-Grashof  case  only,  i.e. 

s  +  l  >  p  +  q.  (6.8) 

Figure  6.3  shows  two  symmetric  configurations  for  the  above  choice  of  parameters. 
We  will  see  that  although  only  two  parameters  Ao  and  A2  are  left  to  be  free,  the  bifurcation 
features  with  respect  to  these  parameters  are  still  informative. 

The  function  now  has  the  following  form: 

=  ^(2A2sm(#32)  +  d2cos(^32)  -  cos(9:i})  +  docos(930) 

+2A 2sm(021 )  -fi  docos(9o\ )  —  dod2co.s(02o)  +  docos(6\o )) 
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Figure  6.3  Symmetric  Configurations 

+  -^-(2.sm(#3o )  +  A2cos(@2o)  —  25m(@io))  ■+■  C  (6-9) 

where  C  is  a  constant  determined  by  d,y  ,  m;  and  the  moments  of  inertia  of  the  bodies, 
I{.  For  i,j  =  0,  1,  2,  3,  6{j  =  9{  —  9j  and  9X  satisfy  the  constraint  equations 

do  -f  cos(@io)  ~p  d2cos(^2o)  d-  cos(8 30)  =  0  (6.10n) 

sin(8\o)  +  d.2sin(62o)  d-  sm(#3o)  =  0.  (6.106) 

In  the  following,  at  symmetric  configuration,  the  variables  will  be  denoted  by 
superscript  “s”  (say,  #j0),  the  formulas  will  be  denoted  by  “|a”  (say,  /(#io)U)-  As 
shown  in  the  example,  the  bifurcation  diagram  of  relative  equilibria  will  be  parameterized 
by  (0,  A0). 

Theorem  6.7:  For  a  floating  four-bar  linkage  with  parameters  shown  in  (6.7),  the 
bifurcation  diagram  of  relative  equilibria  has  the  following  properties: 
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(1)  There  are  two  infinite  branches  in  the  diagram,  one  corresponds  to  0jO  in  the 
leading  form,  another  one  corresponds  to  in  the  lagging  form.  We  refer  to  these  as  the 
symmetric  branches. 

(2)  There  exists  a  constant  Aj£  such  that  no  bifurcation  occurs  on  the  symmetric 
branch  of  leading  form  if  A2  =  A| ;  and,  no  bifurcation  occurs  on  the  symmetric  branch  of 
lagging  form  if  A2  =  —  A£. 

(3)  On  the  symmetric  branch  of  leading  (lagging)  form,  if  A2  <  (~A£) ,  there  exists 

a  constant  c\  (C3)  such  that  the  relative  equilibria  are  stable  for  Ao  <  Ci  (C3),  unstable 
for  Ao  >  Ci  (03),  bifurcated  for  Ao  =  Ci  (03);  one  the  other  hand,  if  A2  >  Aj£  (  —  Aj),  there 
exists  a  constant  c-i  (04)  such  that  the  relative  equilibria  are  unstable  for  Ao  <  C2  (C4), 
stable  for  Ao  >  C2  (c4),  bifurcated  for  Ao  =  ci  (c4). 

(4)  Assume  A2  ^  ±A?j .  Let 


e  =  sgn- 


iei  \\  +  €2  A2  i^3 
e4A2±€5 


(6.11) 


and 


6±  =  sgn(S  1  A2±52) 


(6.12) 


where  ti  and  Si  are  constants  which  are  determined  by  do  and  d'> .  Then,  the  bifurcation 
on  the  symmetric  branch  of  leading  form  will  be  supercritical  pitchfork  if  e+S+  <  0 .  It  will 
be  subcritical  pitchfork  if  e+<5+  >  0.  Similarly,  the  bifurcation  on  the  symmetric  branch 
of  lagging  form  will  be  supercritical  pitchfork  if  e~S~  <0.  It  will  be  subcritical  pitchfork 
if  e~S~  >  0. 

Remarks  6.8: 

(a)  Based  on  the  techniques  of  bifurcation  theory  shown  in  appendix  A,  the  proof  of 
the  above  assertions  is  very  elementary.  However,  two  aspects  have  to  be  considered  before 
the  bifurcation  theory  is  applied.  First,  one  has  to  determine  the  existence  of  a  bifurcation 
and  where  the  branches  are  bifurcated.  Second,  as  one  can  see,  it  is  not  easy  to  write 
down  explicitly  the  function  as  a  function  of  one  variable.  Therefore,  when  applying 
the  techniques  of  bifurcation  theory  which  involves  as  high  as  fourth  order  derivatives  of 
the  function  ,  one  should  consider  the  closure  constraint  equation  (2.2)  simultaneously. 
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(b)  Although  the  proof  of  the  above  assertions  is  elementary,  it  requires  a  large  effort 
in  calculations.  We  used  MACSYMA  to  handle  these  computations.  In  the  following,  we 
only  give  a  sketch  of  the  proof.  | 

Proof  of  Theorem  6.7:  Note  that  the  function  can  be  written  as  a  function  of  relative 
angles  #10,  #20,  and  $30,  which  are  related  through  constraint  equations  (6.10).  From 
(6.10),  we  can  consider  (locally)  #20  and  #30  as  the  functions  of  6 10 .  Again,  from  (6.10) 
one  can  generate  the  quantities  ^  %  I  anj  for  any  positive  integer  i.  Moreover, 

from  Figure  6.3  it  is  easy  to  see  that  at  symmetric  configuration 


t?20  —  7T 


and 


-030- 


(6.13) 


With  above  considerations,  one  can  have  closed  form  expressions  for  '! Jfe- 

^10 

For  instance, 

ffel*  =  icos(°\o)i 

9030  I  -  1 

■9010  Is  _i 


and  fJH*. 


and 


(6.14) 


(6.15) 


and  so  on. 

To  prove  assertion  (1)  in  the  statement  of  the  theorem,  one  needs  to  show  that  at 
symmetric  configuration,  the  equation 

dV{ 


de10  |s 


=  0 


does  not  depend  on  Ao .  Concentrating  on  the  term  involving  Ao  in  and  applying 
(6.13)  and  (6.14),  one  can  show  that  the  first  derivative  of  that  term  with  respect  to  #10 


at  symmetric  configuration  is  zero.  Since  the  rest  of  the  terms  of 


80! 


0  are  still 


functions  of  ,  one  can  see  two  infinite  symmetric  branches  in  the  bifurcation  diagram 
for  two  different  .  Assertion  (1)  is  thus  proved. 

Applying  (6.14)  and  (6.15),  one  can  show  that  the  second  derivative  of  the  function 
at  symmetric  configuration  has  the  form 


Wo 


52n 

dff\Q 


Ao( 


d2  ~  lcosz {e\0)  cos2{6\ 0; 


d,2sin(d*Q ] 


T  A2  ■ 


S 


), 


(6.16) 
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where  II  is  the  summation  of  the  terms  not  involving  Aq  in  .  It  is  obvious  that  when 


x  _  x*  A  2 d2cos2(8s10)  -  d| 

2  '2  co52(^0)«n(^0)’ 


(6.17) 


q 2  ^ 

-r^-L  will  not  depend  on  Ao .  One  can  also  show  that  with  (6.17),  ^  0  under 

assumptions  (6.7)  and  (6.8).  This  means  that  bifurcation  may  not  occur  on  either 
symmetric  branch  of  leading  form  or  symmetric  branch  of  lagging  form.  Note  that  on 
these  different  forms  cos(8l0)  has  the  same  value,  sin(8l0)  has  the  same  absolute  value 
but  different  sign.  (See  Figure  6.3)  Thus,  assertion  (2)  is  proved. 

As  we  have  known  earlier,  the  stability  of  relative  equilibria  depends  on  the  sign  of 
-^r-L.  From  (6.16)  we  see  that  -sp^|s  is  a  linear  function  of  A0 .  Using  the  in  (6.17), 

OVlQ  Ot>10 

the  proof  of  (3)  is  straight  forward. 

To  prove  assertion  (4),  we  apply  the  bifurcation  theory  mentioned  in  appendix  A. 
Let  Aq  denote  c*  in  assertion  (3)  for  some  suitable  i.  One  can  show  that  at  (0io,Ag), 


Moreover 


and 


dVf 

d8io 

d6\0 

d3Ug 

dd'lQd\o 


d2Ve  d3Vc 


dd2o 


d0\ 0  ddxodX 

faa  \*\  ±CiA|  +  e2A2±£3 

(^o,A*)  =  A'1A2±52 


=  0 


(6.18) 


(6.19) 


(6.20) 


for  some  constants  e,  and  6i ,  where  “  +  ”.  corresponds  to  leading  form,  “  —  ”  corresponds 
to  lagging  form.  See  Appendix  B  for  expressions  of  e*  and  Si.  Since  (6.19)  and  (6.20) 
are  not  zero  in  general,  applying  the  lemma  in  appendix  A,  we  can  say  is  strongly 
equivalent  to  the  normal  form  of  pitchfork  bifurcation.  In  addition,  the  type  of  pitchfork 
bifurcation  depends  on  the  sign  of  (6.19)  and  (6.20).  The  assertion  (4)  is  proved.  1 

Remarks  6.9:  (i)  The  condition  of  A2  ^  A^  guarantees  that  (6.20)  and  the  denominator 
of  (6.19)  are  not  zero,  (ii)  In  general  the  bifurcation  changes  from  a  supercritical  one  to 
subcritical  one  at  the  roots  of  numerator  of  (6.19).  1 

Example  6.10;  To  see  how  A2  changes  the  bifurcation  diagram  with  parameter  Ao  we 
give  following  example.  We  will  concentrate  on  the  symmetric  branch  with  respect  to 
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leading  form.  Let  do  =  2  and  d2  —  1 .  Then  e  and  6  have  the  following  form 


e 


+ 


\\  +  6.938A2  —  3.623 
1.873  -  0.051A2 


S+  =  sflrn(1.025  -  0.028A2) 


Then 


e+s+  _  f  -1,  if  -7.426  <  A2  <  0.488; 

(  + 1,  otherwise. 

Note  that  the  region  for  A2  is  an  approximation.  So  we  can  say  that,  when  A2  £ 
(—7.426,  0.488),  the  pitchfork  bifurcation  is  supercritical.  Otherwise,  it  is  subcritical. 
Figure  6.4  shows  this  result.  | 


Before  closing  this  section,  we  would  like  to  make  a  few  additional  remarks. 

(i)  Although  our  discussion  only  concentrated  on  a  structure  of  non-Grashof  type,  a 
version  of  Theorem  6.7  also  holds  for  the  Grashof  case. 

(ii)  Up  to  now  we  have  understood  the  phenomenon  of  bifurcation  on  the  symmetric 
branches.  The  global  analysis  of  the  bifurcations  involves  massive  symbolic  computations. 
However,  as  shown  in  the  Example  6.6  in  this  section  and  other  simulations,  one  can 
numerically  determine  a  global  bifurcation  diagram.  A  large  body  of  such  numerical 
simulations  show  that  the  branches  in  the  bifurcation  diagram  are  connected.  In  other 
words,  for  any  point  in  shape  space,  there  is  a  finite  Ao  for  which  that  shape  determines 
a  relative  equilibrium.  This  is  consistent  with  the  linearity  of  Vf  in  Ao  .  This  property 
provides  a  possibility  to  control  the  attitude  of  a  space  structure  with  a  closed  kinematic 
chain  by  simply  changing  the  position  of  the  centers  of  mass  of  some  bars. 

(iii)  Our  results  in  this  paper  rely  on  some  ideal  conditions,  for  instance,  the 
symmetry  condition  (6.6),  and  the  absence  of  external  and  internal  disturbances.  One 
may  ask  what  will  happen  when  these  conditions  are  violated.  The  answer  to  this  question 
may  relate  to  the  notion  of  universal  unfolding  in  bifurcation  theory.  Numerical  results 
show  that  it  is  possible  to  use  the  unfolding  property  to  control  the  shape  of  the  structure 
near  the  bifurcation  point. 
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7.  GEOMETRIC  PHASES 

Recently,  it  has  become  clear  that  in  a  variety  of  mechanical  contexts,  certain 
(angular)  drift  phenomena  arise  which  can  be  understood  via  geometric  methods.  As 
illustrations,  we  note  the  drift  in  the  orientation  of  a  satellite  caused  by  oscillations  in 
flexible  components  of  the  satellite,  the  slow  change  in  the  plane  of  oscillation  of  a  Foucault 
pendulum  during  a  24-hour  period,  the  ability  of  trampoline  artists  to  reorient  themselves 
in  mid-air  through  limb  and  trunk  movements  etc..  In  the  computation  of  such  drifts, 
one  notes  that  some  or  all  of  the  drift  may  be  attributable  to  geometric  quantities  such 
as  a  path  integral  in  a  space  of  equivalence  classes  (i.e.  shape  space)  of  configurations. 
Hence  the  term  geometric  phase  is  used  in  connection  with  such  drifts.  We  refer  the 
reader  to  (Shapere  &  Wilczek  [1989])  and  (Marsden,  Montgomery  &  Ratiu  [1990])  for 
some  of  the  basic  references  and  the  general  theory.  Of  direct  relevance  to  the  present 
paper  are  the  computation  of  geometric  phases  in  multibody  systems  in  (Krishnaprasad 
[1990];  Krishnaprasad  &  Yang  [1991];  Krishnaprasad,  Yang  &  Dayawansa,  [1991]).  See 
also  (Shapere  [1989])  for  an  example  of  translational  drift  arising  in  swimming  organisms. 
In  the  present  section,  we  will  compute  geometric  phases  for  the  four-bar  linkage  and  show 
special  features  arising  from  the  presence  of  an  closed  kinematic  chain.  We  begin  with  a 
brief  summary  of  the  abstract  framework.  Consider 

(a)  a  simple  mechanical  system  with  symmetry  (Q,  K,  V ,  G) ; 

(b)  a  principal  G -bundle  (Q,Q/G,  G); 

(c)  controls  (forces)  acting  on  (Q,K,  V,  G),  also  leaving  invariant  the  conserved 
momentum  map  J **  :  TQ  — >  Q*  associated  to  the  free  hamiltonian  system  with  energy 
K  +  V,  where  Q*  is  the  dual  of  the  Lie  algebra  Q  of  G.  The  map  J*  is  given  explicitly 
by  the  formula 

where  £  6  Q ,  iG  is  the  Legendre  transform  and  is  the  infinitesimal  generator  (vector 
field  on  Q)  associated  to  £ . 

Define  I?  the  symmetric  bilinear  form  on  Q  by 

I?(£w)  =  ^(£q(<z),7q(<z))- 
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Let  1^  :  Q  — »  Q*  be  the  corresponding  pairing.  Then  we  have  a  vertical-horizontal  splitting 
of  the  tangent  bundle  TQ , 

TQq  =  ( V ert)q  0  (Hor)g 

«,  =  ((I‘)"V)«(?)'+  K  -  ((Ii)-V)o(ff)) 

where  /U  =  J**(u>g).  This  splitting  has  the  equivariance  property  with  respect  to  the  G- 
action  on  Q  and  defines  a  principal  connection.  This  connection,  referred  to  as  mechanical 
connection ,  appears  to  be  originally  due  to  Smale  and  Kummer  (see  (Marsden,  Montgomery 
&  Ratiu  [1990])). 

In  the  concrete  setting  of  planar  many-body  problems  with  rigid  and  flexible 
components,  in  which  G  =  S1  is  abelian,  the  phase  shift  (angle)  associated  to  a  specific 
path  7  in  shape  space  5  =  Q/G  is  simply  the  integral  of  the  connection  form,  over  the 
path.  For  loop  7  it  is  the  holonomy.  In  (Krishnaprasad  [1990]),  a  formula  for  the  phase 
shift  for  a  system  of  N  planar  rigid  laminae  connected  by  (N  —  1)  single  degree  of  freedom 
(N  —  1)  pin  joints  is  derived.  For  this  system,  the  configuration  space  is  TN  parametrized 
by  9  =  (0i,...,0jv)  and  the  shape  space  is  T^”1  parameterized  by  4>  =  (<j>\ , ...,  , 

where  <f>i  =  0,-+1  —  •  Let 


7  =  {(^i(f),  ...<£jv-i(t))|t  E  [0,T],  for  some  T  >  0} 


be  a  closed  curve  in  1 ,  i.e.  <^>;(0)  —  <^,-(T).  Then  absolute  reorientation  or  phase  shift 
of  the  first  body  is  given  by 


A0i 


9(T)  -  0(0)  = 


rT 

f  ~\dt~ 

J 0  e  •  Je 


e  •  J  Md(f> 

e  •  Je 


(7.1) 


where  e  =  (1, ...,  1)T  £  RN  ;  M  £  1)  jg  defined  by 


Mij  =  {  1  if  *  >f  -  1’ 
l  0  otherwise. 

J  is  a  configuration-dependent  N  X  N  matrix,  which  defines  a  riemannian  metric  on  , 

is  invariant  under  the  S'1  action:  (0i,...,0jy)  — >  (0j  +  0,  ...,0/ y  +  0)  for  0  £  F1  (in  other 

words,  it  is  a  matrix  of  functions  of  <j>)\  and,  the  constant  ji  =  e  •  J0  is  angular  momentum 
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of  the  system.  Since  7  is  closed,  A  6{  =  A#i  ,  for  i  =  2 Ar.  Equation  (7.1)  computes 
the  reorientation  of  the  whole  system.  The  first  term  of  (7.1)  is  referred  to  as  dynamic 
phase  which  will  be  zero  if  the  system  is  originally  at  rest;  the  second  term  is  referred  to 
as  geometric  phase  which  will  be  of  interest  to  us. 

We  now  turn  to  systems  of  floating  four-bar  linkages.  As  we  have  seen,  the  riemannian 
metric  given  by  J  in  Section  4  is  also  invariant  under  the  S'1  action.  This  implies  that 
the  formula  (7.1)  is  applicable  to  floating  four-bar  linkages.  We  will  consider  Grashof  and 
non-Grashof  mechanisms  separately. 

Recall  that,  for  Grashof  mechanism,  the  configuration  space  (after  center  of  mass 
reduction)  is  S1  x  S1 ,  which  can  be  parameterized  by  (#0,  #1)  globally.  The  corresponding 
shape  space  is  simply  S1 .  Recall  that  the  angular  momentum  of  the  system  is 

p  =  (l,l)J0  (7.2) 

where  0  —  (<?o,0i).  Since  J  is  a  matrix  of  functions  of  #10  =  <j> ,  one  can  derive  the  same 
equation  as  (7.1)  with  N  =  2.  By  appealing  to  Proposition  4.3,  one  gets 

A#0  =  /  A-cfc  -  f  —  0  (7.3) 

Jo  e  •  Je  J1  e  •  Je 

where  e  =  (l,l,l,l)T,  J  has  been  shown  in  Section  4  and  M  is  the  second  column  of 
(I,Q,T)T  matrix,  i.e. 


jdfT  _  (  r>  1  _  lisin(93  —  8 1)  hsiiUJh^JL l)  f 
\  l2sin(03  —  $2)  l3sin(03  —  02)  ' 

Recall  that  in  Grashof  case,  the  configuration  space  has  two  connected  components, 
corresponding  to  leading  form  and  lagging  form,  respectively.  On  these  components,  J  has 
different  symbolic  expressions.  Consequently,  the  phase  shift  generated  by  the  same  closed 
curve  in  shape  space  will  be  different.  We  illustrate  this  point  by  the  following  example. 

Example  7.1: 

Let  ji  =  0.  We  choose  the  parameters  as  follows. 

mo  =  1,  m\  —  1,  m2  =  1,  m3  =  1; 
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io  1,  I\  —  1,  I2  —  1,  J3  —  1; 


d03  =(-1.5,0),  d01  =(1.5,0),  d10  =  (-0.5,0),  d12  =  (0.5,0), 

d21  =  (-1.5,0),  d23  =  (1.5,0),  d32  =  (-2,  -1),  d30  =  (2,  -1). 


Figure  7.1 

Figure  /.I  shows,  for  both  leading  and  lagging  form,  how  the  orientation  of  the  0-th 
bar  varies  when  the  1-st  bar  revolves  about  the  0-th  bar.  One  can  see  that  when  the 
linkage  returns  to  its  original  shape  (in  this  example,  the  0-th  bar  and  1-st  bar  are  folded 
together),  the  0-th  bar  has  been  re-oriented  by  A$o  which  is  less  than  zero.  Its  value 
for  leading  form  is  different  from  the  one  for  lagging  form  because  the  existence  of  the 
mass-offset  in  the  3-rd  bar.  Here,  we  see  again  in  Grashof  mechanisms,  two  forms  should 
be  considered  separately.  g 

Now  we  consider  non-Grashof  mechanisms.  Recall  that  in  this  case,  the  configuration 
space  is  still  51  x  S'  and  the  symmetry  group  is  F1  .  However,  the  physical  angles  6q 
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and  6\  cannot  parameterize  the  configuration  space  globally.  Nevertheless,  the  angular 
momentum  of  the  system  can  still  be  expressed  as  (7.2)  locally.  More  precisely,  every 
element  in  the  matrix  J  is  a  branched  function  (branched  by  corresponding  leading  and 
lagging  form)  and  not  defined  at  the  angle  limits  of  6\q  .  Since  these  limits  are  only  two 
isolated  points  in  shape  space  S'1 ,  the  integration  (7.3)  can  be  performed.  Let  7  be  the 
closed  loop  on  Sl  given  by 

7  =  {6\o(.t)\t  e  Ai}  U  {0io(*)|<  G  A2}  =  7 lead  u  7 lag 


where 


Ai  =  {t  e  [0, T]\sin(dz2)  <  0(  leading  form}) 

A 2  -  {t  G  [O,T]|sm(032)  >  0(  lagging  form}) 
and  Ai  U  A2  =  [0,  T] .  Then  the  phase  reconstruction  (net  reorientation)  is  given  by 


A0n  = 


[  -J^dt+  [  -JL-dt-  [ 

j  A  1  ^  '  J  ^  J ^  ‘  d  Q  J  ~)'le  ad 


e- J M 
e  •  Je 


dd 


10 


flag 


e-JM 
e  •  Je 


d9io 


(7.4) 


In  the  following  we  show  an  example  of  the  non-Grashof  case. 

Example  7.2: 

Again,  let  ji  —  0.  The  parameters  are  given  as  follows. 


m0  =  5,  m  i  =  l, 

m2  =  1, 

m3  =  1; 

Iq  =5,  I\  =  1, 

/2  =  1, 

h  =  1; 

do3  =  ( — 2, 0),  d0i  =  (2,0), 

dio  =  ( 

-1.5,0), 

di2  =  (1-5,0 

d2i  =  (—1.5,  A),  d23  =  (1.5,  A), 

d.32  =  1 

(-1.5,0), 

d30  =(1.5,0). 

It  should  be  noted  that  this  system  has  symmetric  configurations  for  any  A .  The  scheme 
to  change  the  shape  of  the  system  is  given  by  the  following  steps. 

1.  Starting  at  a  symmetry  configuration  in  leading  form; 

2.  Increasing  #10  —  tt  until  it  reachs  the  dead  point; 

3.  Decreasing  9\ 0  —  tx  in  lagging  form  until  it  reaches  the  second  dead  point; 

4.  Increasing  01O  —  tt  in  leading  form  until  it  reaches  the  original  symmetric  configura¬ 
tion. 
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Figure  7.2  shows  the  change  of  9q  versus  the  change  of  the  shape  of  the  linkage  by  the 
above  scheme  for  different  A .  It  is  interesting  to  observe  that  the  value  for  re-orientation  of 
the  0-th  bar  can  be  positive,  negative  or  zero.  One  can  show  that  when  a  floating  four-bar 
linkage  is  at  rest  at  its  symmetric  configuration,  varying  the  position  of  center  of  mass  of 
the  0-th  or  the  2-nd  bar  so  that  the  condition  (6.6)  is  unchanged  will  not  change  the  shape 
and  the  orientation  of  the  system.  Then,  this  example  suggests  that  the  system  can  be 
reoriented  to  arbitrary  direction  by  allowing  the  center  of  mass  of  the  0-th  or  2-nd  bar  to 
be  movable  so  that  the  condition  (6.6)  is  kept  satisfied.  | 
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8.  CONCLUSION 

The  kinematics  and  dynamics  of  planar,  floating  four-bar  linkages  have  been  investi¬ 
gated  in  this  paper.  We  have  observed  that  the  topology  of  the  configuration  space  for  such 
linkages  is  related  to  the  classical  Grashof  criterion.  We  have  shown  that,  depending  on 
the  ratios  of  the  lengths  of  the  linkages,  the  configuration  space  is  the  union  of  two  copies 
of  R?  x  S1  x  51  ’s  (Grashof  mechanism),  or  a  single  copies  of  R?  x  51  x  Sl  (non-Grashof 
mechanism),  or  R2  x  S1  x  “ figure  8”  (change-point  mechanism)  which  is  not  a  smooth 
manifold.  We  have  derived  a  closed  form  of  the  lagrangian  for  the  system  after  account¬ 
ing  for  the  loop  closure  condition  (potential  energy  has  been  ignored).  The  expression  for 
the  lagrangian  reveals  that  systems  of  floating  four-bar  linkages  possess  S1  symmetry  after 
center  of  mass  reduction.  This  property  made  it  possible  for  us  to  understand  the  dynamics 
of  the  system  via  the  theory  of  symplectic  and  Poisson  reduction.  Although  the  geometric 
structure  of  this  system  is  parallel  to  that  of  the  planar  two-body  problem  (Sreenath  et 
al.  [1988]),  because  of  the  presence  of  kinematic  loop  constraints  and  greater  freedom  in 
the  choice  of  structural  parameters,  floating  four-bar  linkages  have  more  complicated  and 
interesting  dynamics.  We  have  demonstrated  this  point  by  studying  the  relative  equilibria 
for  those  linkages  which  admit  symmetric  configurations.  As  we  have  shown,  a  symmetric 
configuration  is  itself  a  relative  equilibrium.  If  we  vary  certain  kinematic  parameters  which 
preserve  the  symmetry,  a  symmetric  relative  equilibrium  may  be  bifurcated.  The  type  of 
bifurcations  can  be  either  supercritical  or  subcritical  pitchfork.  The  stability  of  the  relative 
equilibria  at  symmetric  configurations  was  also  investigated.  The  examples  have  shown 
that  a  floating  four-bar  linkage  can  have  as  many  as  ten  relative  equilibria.  Finally,  we 
investigated  the  re-orientation  problem.  We  have  shown  that  the  formula  for  calculating 
the  geometric  phase  for  an  open  chain  system  can  be  extended  to  floating  four-bar  linkages. 
Through  an  example,  we  showed  that  for  non-Grashof  mechanisms  which  admit  symmetric 
shape,  the  geometric  phase  shift  may  be  generated  in  arbitrary  direction. 

As  we  have  mentioned  in  the  first  section,  a  floating  four-bar  linkage  is  the  simplest, 
but  non-trivial  multibody  systems  with  closed  loop.  Fully  understanding  the  kinematics 
and  dynamics  of  these  kinds  of  systems  will  certainly  be  useful  in  investigating  more 
complicated  systems.  One  can  expect  that  the  kinematic  and  dynamic  properties  of  floating 
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four-bar  linkages  discovered  in  this  paper  will  also  show  up  in  more  complex  systems.  This 
will  be  pursued  in  future  work. 


45 


Appendix  A. 

In  this  appendix,  we  review  some  basic  concepts  and  results  of  bifurcation  theory. 
The  standard  reference  is  (Golubitsky  &  Schaeffer  [1985]). 

Consider  a  single  scalar  equation 

<?(x,A)  =  0  (-4.1) 


The  bifurcation  theory  studies  how  the  solutions  x  of  this  equation  change  with  the 
parameter  A ;  or,  more  precisely,  what  type  of  bifurcation  occurs  with  parameter  A . 
Without  loss  of  generality,  one  can  assume  <?(0, 0)  =  0 .  Moreover,  we  assume  g  :  Rx  R  — >  R 
is  smooth.  This  is  one  of  the  standard  local  (static)  bifurcation  problems  with  one  state 
variable,  called  the  recognition  ■problem.  As  with  many  bifurcation  problems,  this  problem 
can  be  solved  successfully  through  singularity  theory,  in  which  the  related  issue  is  called 
finite  determinancy.  Let  2  =  (a:,  A).  Near  origin,  the  function  g  can  be  written  as 


9W=  £  h{Yzrg{0)Z°+  £  “"W2”’  (A2) 

|a|<fc+l  ‘  |a|=fc+l 

for  some  smooth  functions  aQ  defined  in  a  neighborhood  of  the  origin.  Here  we  used  the 
conventions  with  multi-indices: 


|a|  -  op  +  a2,  a!  =  (ai)!(a2)!, 

o  O'-) 

__  a-1  \  a2  ( _\ot  _  (  °  \ct\(  °  \ot2 

’  ^dz’  ' dx >  ' 

A  key  question  is  what  terms  in  (A. 2)  can  be  ignored  such  that  the  values  of  coefficients  of 
remaining  terms  can  be  used  to  determine  the  qualitative  behavior  of  the  original  equation 
(A.l),  for  example,  the  variation  in  number  of  solutions.  Singularity  theory  solves  this 
problem  by  finding  a  suitable  change  of  coordinates  such  that  function  g  is  equivalent 
to  a  standard  model  h ,  called  normal  form.  A  precise  definition  is  given  as  follows  (see 
(Golubitsky  &  Schaeffer  [1985])). 

Definition  A.l:  Two  smooth  mappings  g,h  :  Rx  R  — »  R  defined  near  the  origin  are 
equivalent  if  there  exist  a  local  diffeomorphism  of  R2 ,  (x,  A)  (A(x,  A),  A(A))  at  the 
origin  and  a  nonzero  function  5(x,  A),  such  that 

g(x,  A)  =  S(x,  A )h(X(x,  A),  A(A))  (A3) 
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where  ^(0,0)  >  0  and  Ar(0)  >0.  If  A  =  A,  g  and  h  are  strongly  equivalent.  y 

From  this  definition  we  see  that,  since  S(: x,  A)  is  nonzero,  the  solution  of  g{ x,  A)  =  0 
and  h(X,  A)  =  0  are  the  same  in  the  sense  of  diffeomorphism.  From  this  point  of  view, 
by  means  of  singularity  theory  one  can  show  why  and  what  the  high-order  terms  in  (A. 2) 
do  not  effect  the  qualitative  behavior  of  equation  g(x,  A)  =  0.  It  should  be  noticed  that 
although  this  method  does  not  tell  us  how  to  derive  an  appropriate  normal  form  h ,  for 
most  physical  problems,  such  as  the  one  considered  in  this  paper,  it  is  not  hard  to  pick  up 
some  of  the  candidates  from  a  large  number  of  known  simple  polynomials  of  x  and  A ,  or 
the  model  of  normal  forms  which  have  standard  bifurcation  diagrams.  This  is  in  essence 
the  spirit  of  application  of  singularity  theory  to  a  physical  problem. 

Without  considering  detailed  issues  of  singularity  theory  which  are  applicable  to 
bifurcation  problems,  we  directly  give  the  following  result  which  will  be  used  in  the  next 
section.  For  details  see  (Golubitsky  &  Schaeffer  [1985]),  chapter  2.  First  we  need  the 
concept  of  germs. 

Definition  A. 2:  Two  smooth  functions  defined  near  the  origin  are  equivalent  as  germs  if 
there  is  some  neighborhood  of  the  origin  on  which  they  coincide.  Let  £x,\  denote  the  set 
of  equivalence  classes  of  such  functions.  The  elements  in  £X:\  are  called  germs.  | 

Lemma  A. 3:  A  germ  g  G  £x.\  is  strongly  equivalent  to 

exk  +  6\x  (A. 4) 


for  k  >  2  if  and  only  if  at  x  =  A  =  0 


d 

,J  =  TJ 


and 


,  d  c  d  d 

£  =  s9n(-)9,  S=sSn-^g. 


(A.  5a) 


{A.  56) 


Remark  A. 4:  When  k  —  3  the  normal  form  (A. 4)  provides  a  pitchfork  bifurcation.  From 
this  lemma,  so  is  g  if  (A. 5)  holds.  It  is  easy  to  show  that  if  e8  >  0 ,  the  pitchfork  bifurcation 
is  subcritical;  if  eS  <  0 ,  the  pitchfork  bifurcation  is  supercritical.  flj 
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Appendix  B 

The  expression  of  and  Si  in  the  theorem  6.7  are  of  the  following  form. 


ei 

=  — 12do(d2 

—  do)2(d2  —  dod2  —  d2d2  —  2d2 

+  d30- 

4do); 

^2 

=  48d2d2((d; 

2  —  do)2  —  l)(d2  —  do)\/ 4  —  (d( 

)  —  d2)2 

=  48d2d2(d2 

—  do)2(d2  —  dod2  —  dgd2  —  4d2 

;  +  do  — 

2do); 

=  d2((d2  —  do)2  —  4)(d2  —  do)2  yj  A  —  (do  — 

d2)2; 

^5 

=  2^((d2  - 

d0)2  —  4)((d2  —  d0)3  —  4d2). 

A  =  4d2 \A  —  (do  —  d2  )2 . 


Then  i 

51  =  -^(d2  -  d0)V4-(d0-d2)2; 

2 

52  —  — — d2(d2  +  3dod2  —  3cf q c?2  4~  4d2  T  dj)). 
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